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Abstract 

Wc report on the results of detailed spectral studies of the ultra- luminous X-ray source (ULX), M81 
X-9 (Holmberg IX X-l), made with XMM-Newton on 2001 April 22 and with ASCA on 1999 April 6. On 
both occasions, the source showed an unabsorbed 0.5-10 keV luminosity of ~ 2 x 10 40 erg s _1 (assuming 
a distance of 3.4 Mpc) and a soft spectrum apparently represented by a multi-color disk model with an 
innermost disk temperature of 1.3—1.5 keV. Adding a power-law model further improved the fit. However, 
as previously reported, the high luminosity cannot be reconciled with the high disk temperature within a 
framework of the standard accretion disk radiating at a sub-Eddington luminosity. Therefore, we modified 
the multi-color disk model, and allowed the local disk temperature to scale as oc r~ p on the distance r 
from the black hole, with p being a free parameter. We then found that the XMM-Newton and the ASCA 
spectra can be both reproduced successfully with p ~ 0.6 and the innermost disk temperature of 1.4-1.8 
keV. These flatter temperature profiles suggest deviation from the standard Shakura-Sunyaev disk, and 
are consistent with predictions of a slim disk model. 

Key words: accretion, accretion disks — black hole physics — X-rays: individual (M81 X-9, Holmberg 
IX X-l) — X-rays: stars 



1. Introduction 

Nearby spiral galaxies often host ultra-luminous com- 
pact X-ray sources (hereafter ULXs; Makishima et 
al. 2000), with X-ray luminosities reaching 3 x 10 39 - 
10 40 erg s~ 4 . Although they require black hole masses of 
20-100 Mq or more by considering the Eddington limit, 
Le, it is not obvious whether such massive stellar black 
holes can be created as an end point of evolution of mas- 
sive stars. Thus, there are still extensive arguments on the 
nature of the ULXs. Thanks to ASCA, a great leap has 
been achieved in our understanding of 1-10 keV spectra 
of a dozen ULXs. As reported by many authors (Okada et 
al. 1998; Makishima et al. 2000, and references therein), 
spectra of a majority of the ASCA sample have been 
reproduced by so-called multi-color disk (MCD) model 
(diskbb in xspec 1 ; Mitsuda et al. 1984; Makishima et 
al. 1986), while others are described with a power-law 
model. The MCD model approximates emission from 
a geometrically thin, radiation dominated standard ac- 
cretion disk (Shakura & Sunyacv 1973), and success- 
fully describes the dominant soft spectral component of 
Galactic/Magellanic black hole binaries (hereafter BHB) 

1 http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/ 



in the high/soft state. In addition, spectral transitions 
between the disk dominant type and the power-law domi- 
nant type states have been observed in several ULXs (e.g., 
Kubota, Done, Makishima 2002; La Parola et al. 2001). 
These spectral characteristics of ULXs have thereby rein- 
forced their candidacy as accreting black holes. 

In spite of the successful application of the MCD model 
to the ULX spectra, there still remain three self incon- 
sistencies. One is known as "too high Ti n " problem that 
the values of the innermost disk temperature, Tj n , ob- 
tained from the disk-dominant ULXs (typically 1-2 keV; 
Makishima et al. 2000) are significantly higher than those 
of BHBs (0.5-1.2 keV), even though higher-mass black 
holes should have lower values of T; n as T; n oc Af -1 / 4 under 
the standard-disk picture (e.g., Makishima et al. 2000). 
The second one is known as "variable n n " problem: the 
inner disk radius, r- mi of ULXs are reported to vary as 
r m ocTj" 1 (Mizuno, Kubota, Makishima 2001), in contra- 
diction to the case of high-state BHBs in which r m is gen- 
erally constant at a value consistent with the radius of last 
stable orbit around the central black hole (e.g., Makishima 
et al. 1986; Tanaka & Lewin 1995; Ebisawa et al. 1993; 
1994). The last problem is that the ULXs luminosity at 
which the spectral transition occurs, ~ 10 39 erg s _1 , is too 
high to be considered as the classical high-low state tran- 
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sitions of BHBs which generally occurs at several percent 
of Le 

The "too high Tj n " problem might be solved within the 
framework of standard accretion disks, by invoking ex- 
treme assumptions. For example, a rapid black hole ro- 
tation could increase T- ln (Makishima et al. 2000) but the 
effect may not be sufficient to explain the ULX spectra 
(Ebisawa et al. 2003). Alternatively, the values of T- m 
could increase if the color-to-cffectivc correction factor k 
took a very high values as 3-5, for some unknown reasons, 
but this would contradict the canonical value of k=1. 5-2.0 
suggested theoretically (Shimura & Takahara 1995; Davis 
et al. 2005) and observationally (Makishima et al. 2000). 
Furthermore neither of these "solutions" may solve the 
second and third problems. 

Considering these difficulties, accretion flows in ULXs 
may be significantly deviated from the description of stan- 
dard accretion disks, even though the emergent spectra 
look rather alike. Actually, when the mass accretion rate, 
M, increases and the luminosity becomes closer to Le, the 
disk structure is expected to change from the standard 
disk to a slim disk (e.g., Abramowicz et al. 1988; Watarai 
et al. 2000). The slim disk is predicted to have a signifi- 
cantly higher inner temperature than the standard accre- 
tion disk, and its apparent r- ln (when fitted with the MCD 
model) is expected to vary roughly as oc T^ 1 (Watarai 
et al. 2000) just as is observed. These explain away 
the first and second problems, respectively (e.g., Mizuno, 
Kubota, Makishima 200f; Watarai, Mizuno, Mincshigc 
200f). Further, invoking the slim disk implies that the 
radiation luminosity is close to Le, thus solving the third 
problem. Therefore, the slim disk scenario can potentially 
solve the three problems associated with the interpreta- 
tion of disk dominated ULXs. 

While the radiation-dominated standard disk is char- 
acterized by a local temperature varying as T{r) oc r~ 3 / 4 
(r being the distance from the central black hole), the 
temperature profile of the slim disk is predicted to be 
flatter than this as T(r) oc r^ 5 . The different temper- 
ature profiles are expected, in turn, to give subtle dif- 
ferences to the overall disk spectra. This effect may be 
described by so-called p-free disk model (a local model 
for xspec 2 ; Mineshige et al. f994 ; Hirano et al. f995 
; Kubota & Makishima 2004), it assumes that the spec- 
trum is composed of multi-temperature blackbody emis- 
sion, where the local temperature T(r) at a radius r is 
given by T{r) = T m (r / n n ) _?> , with p being a positive 
free parameter. Using the ASCA and RXTE data of 
Galactic/Magellanic BHBs, it has been confirmed that 
p takes a value consistent with the standard value, 3/4, 
when the objects are considered to reside in the standard- 
disk state with constant values of Hn (Kubota et al. 2005), 
and that p becomes smaller than 3/4 in a supposedly slim 
disk state (Kubota & Makishima 2004; Abe et al. 2005). 

In the present paper, we apply the p-free disk model 
analysis to X-ray spectra of M8I X-9 in the disk dom- 
inant state, obtained with XMM-Newton on 2000 April 

2 http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspcc/modcls/diskpl 



22 and with ASCA on 1999 April 6, in search for the ex- 
pected spectral deviation from the standard disk. At a 
distance of 3.4 Mpc (Georgicv et al. 1991), this source 
is an extensively studied ULX associated with the dwarf 
galaxy Holmberg IX, which is a close companion to M8I. 
Through previous snapshot observations with ASCA in 
1993-1999, its unabsorbed 0.5-10 keV luminosity was 
found in (6-20) xlO 39 erg s _1 assuming an isotropic emis- 
sion. Until 1998, the source kept a relatively low luminos- 
ity of ~ 6 x 10 39 erg s _1 , with a single power-law spectrum 
of photon index 1.4 (Ezoc, Iyomoto, Makishima 2001). 
In the 1999 observation, the source was reported to have 
made a spectral transition to a disk-dominant state with 
Ti n = 1-24 keV, together with an increased luminosity of 
- 1.8 x 10 40 erg s" 1 (La Parola et al. 2001). Good XMM- 
Newton observations of M81 X-9 were performed three 
times by 2002. In 2001, its luminosity and spectral shape 
were almost the same as those obtained in the 1999 ASCA 
observation (Sugiho 2003). In the 2002 observations, the 
source became slightly fainter with an unabsorbed 0.3- 
10 keV luminosity of L = 1.1-1.3 x 10 40 erg s _1 , and its 
spectrum returned to the power-law dominant state with 
a hint of cool MCD component (Miller, Fabian, Miller 
2004). 

Out of these multiple datasets, here we select two par- 
ticular ones for our p-free study, namely the 1999 ASCA 
data and the 2001 XMM-Newton data, because the object 
was clearly in the disk dominant state on these occasions. 
These data sets are expected to give a significantly smaller 
values of p than that of standard disk, 3/4, if the slim disk 
is realized. These are the same data as analyzed already 
by La Parola et al. (2001) and Sugiho (2003), respectively. 

2. Observations and Data Reduction 

The present XMM-Newton observation of M81 X-9 was 
performed for 127.9 ks on 2001 April 22 (observation ID 
of 0111800101), with the three detectors of the European 
Photon Imaging Camera (EPIC). The two MOS cameras 
were operated in full frame mode, and the PN camera with 
the medium filter was operated in small window mode. We 
retrieved the data from the XMM-Newton public archive 3 , 
and used only MOS-1 data for the spectral study, since the 
source fell in the chip boundary of MOS-2 and outside the 
field of view (FOV) of the PN camera. The event file was 
reprocessed and filtered with the XMM-Newton Science 
Analysis Software (SAS v6.0.0). Utilizing light curves of 
source-free regions from the MOS-1 data, we rejected in- 
tervals with high background levels, and obtained a net 
exposure of 79.7 ks for the MOS-1 camera. 

Figure 1 shows a MOS-1 image in the range of 0.5-10 
keV. The source is thus located at 12.'5 off the center of 
the MOS-1 FOV. Using a SAS task xmmselect (v6.0.1), 
the MOS-1 events with patterns to 12 were extracted 
from a circular region of 50" radius around the position of 
M81 X-9, (a, 6) = (09 h 57 m 53 s , +69°03'47"). Background 
events were accumulated over a source-free region with a 

.fitmlhttp:/ /xmm. vilspa.esa.es/ 
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I. '7 radius, at a position (a, 6) = (09 h 57 m 43 s ,+68°59'24"). 
The on-source and the background regions are indicated 
in figure 1. Figure 2 shows the 0.5-10 keV light curve of 
M81 X-9. There is no significant intensity variations be- 
yond the Poisson noise, with a typical upper limit of 10%. 
The average count rate, 0.424 ±0.004 cts s _1 , roughly cor- 
responds to a 0.5-10 keV flux of 1.4 x 10~ n erg s _1 cm~ 2 . 

Since M81 X-9 is located close to the edge of MOS-1 
FOV, its spectral analysis might be subject to calibration 
uncertainties. Therefore, we analyzed three data sets of a 
calibration target, the supernova remnant G2 1.5-09 which 
shows a power-law spectrum, one acquired near the FOV 
center and the others 1 1 ' 1 2 and 13^02 off. As detailed in 
Appendix, we have confirmed that the spectra of such off- 
axis sources can be studied with a high reliability. Though 
the off-axis data may yield slightly harder spectrum than 
the on-axis ones, the difference in the photon index is less 
than ~ 0.1 (see figure6 in Appendix). 

To compare with the XMM-Ncwton result on the disk 
dominant state, we also reanalyzed the ASCA SIS (Solid 
state Imaging Spectrometer) data of M81 X-9 obtained 
on 1999 April 6. We used the screened data set pro- 
vided by the DARTS science archive 4 with a net expo- 
sure of 33 ks. The source spectra were accumulated over 
a circular region centered on the source, with a radius 
of 2 '9, and the background spectra were created using 
blank sky data. The criteria of the data screening are 
as follows: a) the object be at least 10° above the night 
Earth's limb or at least 20° above the bright Earth's limb; 
c) the cutoff rigidity (COR) of cosmic rays be greater than 
6 GeV c^ 1 ; and d) the time after day night transition be 
greater than 100 s. The 0.5-10 keV count rate of M81 X-9 
was 0.326±0.003 and 0.272±0.003 cts s" 1 with SIS0 and 
SIS1, respectively. The implied 0.5-10 keV flux, approxi- 
mately 1.3 x 10~ n erg s _1 cm~ 2 , is close (~ 93 %) to that 
of the 2001 XMM-Newton data. 

3. Data Analysis and Results 

3.1. Analysis of the XMM-Newton data with canonical 
models 

A left panel of figure 3 shows a MOS-1 spectrum of 
M81 X-9 in the range of 0.5-10 keV. Using a SAS task 
epatplot, the MOS-1 data were confirmed to be free 
from photon pile-up effects. Following the standard mod- 
eling of ULXs (e.g., Makishima et al. 2000), we first 
examined the spectrum with a single power-law model 
and the MCD model, and obtained the fitting results as 
shown in Table 1. Here, the response matrices were cre- 
ated by rmfgen and ARFGEN (v6.0.1), and the spectral 
fitting was performed with the software package XSPEC 

II. 2.0. The model was subjected to photoelectric ab- 
sorption by a hydrogen column which is left free to vary. 
The single power-law fit was far from acceptable, with 
X 2 /(d.o.f) = 574.1/335. The single MCD fit was much 
more successful, and yielded the parameters which are 
very similar to the ASCA results in 1999 reported by La 




Fig. 1. A 0.5-10 keV unsmoothed MOS-1 image of the M81 
X-9 field obtained on 2001 April 22. The source region and 
the background region are shown with a solid circle and a 
dotted circle, respectively. The background is included. 




2x10* 3x10* 4x10* 
TIME s 



5x10* 6x10* 



Fig. 2. The 0.5-10 keV MOS-1 light curve of M81 X-9, plot- 
ted with a time bin of 1 ks. A solid line indicates the best fit 
constant. 
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Parola et al. (2001). However, the fit still remained unac- 
ceptable with x 2 /( d -°- f ) = 427.4/335. 

We hence fit the data with an MCD plus power-law 
model, as is done for galactic BHBs in the high state. This 
two-component model has successfully reproduced the ob- 
served spectrum with x 2 1 (d.o.f) = 359.4/333. As shown in 
Table 1, the best-fit model implies an absorbed 0.5-10 keV 
flux of 1.4 x 10~ n erg s _1 cm~ 2 , in which the MCD com- 
ponent accounts for ~ 60%. The unabsorbed luminosity in 
the same band is obtained as 2.3 x 10 40 erg s _1 , assuming 
an isotropic emission. The MCD parameters, Tj n and n n , 
have been measured to be 1.45 keV and (160 ± 12) • a km 
with a = -\/cos60 / cosz (i being the inclination angle of 
the disk) , respectively. The true inner radius of the disk 
is calculated as R{ n = (190 ± 15) • a km, by using a rela- 
tion of R[ n = (ft 2 ri n , when £ = 0.41 is a correction factor 
for the inner boundary condition (Kubota ct al. 1998) 
and k — 1.7 is a spectral hardening factor (Shimura & 
Takahara 1995). The disk bolometric luminosity is esti- 
mated as Lboi = 1-5 x 10 40 • a 2 erg s _1 through a relation 
of Lboi = 47TCTr? n J; 4 (Makishima et al. 1986), with a the 
Stefan-Boltzman constant. When the 0.5-10 keV unab- 
sorbed power-law luminosity, 8.7 x 10 39 erg s _1 , is added 
to the disk bolometric luminosity, the total luminosity is 
roughly estimated as 2.4 x 10 40 erg s" 1 for a = 1 (i.e., 
i = 60°). 

The observed X-ray luminosity requires a black hole 
mass of > 150 M Q if the Eddington limit is satisfied. In 
contrast, the observed disk temperature of 1.45 keV is too 
high for the emission to be coming from a standard disk 
around such a massive black hole. In other words, the 
disk inner radius of 190 km is too small to be identified 
with the last stable orbit of a non-spinning black hole with 
> 150 M . If the obtained value of i?j n were required to 
coincide with 6R g , then the Schwarzshild radius, 2i? g , of 
this black hole would be only ~ 64 km, which is equivalent 
to a black hole mass of ~ 20 M . Thus, the 2001 XMM- 
Newton observation of M81 X-9 clearly reproduces the 
"too high Tin" problem described in § 1. 

3.2. Analysis of the ASCA data in 1999 

The flux and the MCD parameters obtained from the 
2001 XMM-Newton data (§ 3.1) are very similar to those 
reported by La Parola et al. (2001), who analyzed the 
ASCA SIS data obtained on 1999 April 6 with the sin- 
gle power-law model and the single MCD model as in- 
troduced in § 3.1. In order to compare the ASCA and 
XMM-Newton spectra in a unified manner, we reanalyzed 
the same 0.5-10 keV ASCA SIS data with the same two- 
component model. 

A right panel of figure 3 shows the ASCA SIS (SIS0 and 
SIS1) spectra with the best fit models. The results of this 
analysis are also shown in Table 1. The obtained best-fit 
parameters of the single power-law fit and the single MCD 
fit are consistent with those reported by La Parola et al. 
(2001). Although the single MCD model gave an accept- 
able fit with x 2 1 (d.o.f) = 318/301, addition of a power-law 
component improved the fit significantly in terms of an F- 
test, with F(2,301) = 20.8 which corresponds to 99.99 % 



significance. This two-component fit implies a 0.5-10 keV 
flux of 1.3 x 10 _11 crg s~ 1 cm -2 with the MCD component 
contributing ~ 70%, and an unabsorbed 0.5-10 keV lumi- 
nosity of 2.4 x 10 40 erg s _1 . The values of T- ln and r- m have 
been estimated as 1.28±g;^| keV and (220 ±50) -a km, 
respectively (Table 1). The spectral shape and the lumi- 
nosity are thus almost the same, within errors, between 
the two observations with the different instruments. 

3. 3. Temperature profile of the optically thick disk 

In addition to the "too high T; n " problem, we en- 
countered another problem in the spectral fit of M81 X- 
9. Although the two-component (MCD plus power-law) 
model apparently reproduces the spectra from the ASCA 
and the XMM-Newton observations, the power-law com- 
ponent exceeds the MCD component at energies below ~ 
1 keV, as shown in figure 3. Such a fit would be rather un- 
physical, since the power-law tail is generally considered as 
a consequence of inverse- Compton scattering of the MCD 
photons by high-energy electrons created around the ac- 
cretion disk. Therefore, the hard-tail flux should not dom- 
inate that of the seed MCD, at least in the energy band 
below 2.8 fcTj n which is the peak energy of a blackbody 
emission of temperature T; n . The apparent low-energy 
dominance of the power-law component suggests a sub- 
tle deviation of the accretion disk configuration from the 
standard one, in such a way that the softest end of the 
radially-integrated spectrum is more enhanced than the 
MCD model. 

As described in §1, such a spectral change may arise if, 
e.g., the radiative efficiency of the inner disk region be- 
comes reduced and the radial temperature gradient flat- 
tens; such effects can be quantified using the p-free disk 
model (Kubota & Makishima 2004). We hence fitted the 
MOS-1 and SIS data with the p-free disk model, and ob- 
tained the results presented in Table 1 and figure 3. In 
both the XMM-Newton and ASCA spectra, the p-free 
disk model is similarly, or even more, successful than the 
MCD+power-law fit. The values of p have been found as 
0.60±0.02 and 0.61+g;g| with the Newton and ASCA data, 
respectively, in a very good mutual agreement. They are 
significantly smaller than the standard value of p = 3/4. 
The absorption column becomes somewhat higher than 
that implied by the single MCD fit. To reinforce the re- 
sult, we simulated the spectrum by considering the possi- 
ble systematic difference in the spectral index of Ar w 0.1 
for the Newton off-axis sources (see figure 6 in Appendix). 
This difference was found to make the best fit value of p 
slightly smaller as p = 0.58 ± 0.02, but it is within the 
statistic uncertainty. 



No. 



Accretion disk in M81 X-9 



Table 1. Spectral fitting parameters and derived quantities of M81 X-9. 



Model JVg r or pt T in (keV) r in (km)S M M cd (A/ ) ss -Fx" ^boi # xVd.o.f. 

2001 April 22 (XMM-Ncwton) 

574.1/335 
427.4/335 
359.4/333 

357.3/334 



power-law 

MCD 

MCD 

+powcr-law 
p-free disk 



n yc+0.08 
u - '°-0.07 
i 7 +0.63 
/ - 1 ' -0.42 



1 7n+ 019 

1 - ' u -0.20 



1.9610.03 



9 90+0. 58 
Z - ZO -0.35 

0.60+0.02 



1.35 
1.45 



0.02 

0.03 
+0.10 
0.14 



1 79+0.11 
1 - > Z -0.10 



220 ± 4 
160 ± 12 

100 ± 20 



30 
22 



14.9 

13.0 21.0 

8.4 14.6 
5.1 

13.5 28.4" 



1999 April 06 (ASCA) 



power-law 

MCD 

MCD 

+power-law 
p-frec disk 



4-68^ 
1.71+0.14 



2.1610.04 



3.18 



+2.76 
-0.94 



1.2510.03 
1.28 



2.68 



+0.37 

-0.38 



2.35 
0.61 



1.95 

0.65 
+0.04 
0.03 



+0.12 
-0.14 



260 
220 



13 
50 



35 
29 



1.47 



+0.12 
-0.10 



140 ± 30 



13.9 
12.4 
8.5 
4.3 

12.6 28.3" 



21.8 
16.1 



453.3/301 
318.5/301 
297.8/299 

298.3/300 



* All quoted uncertainties are at the 90 % confidence limit. 

t Hydrogen column density for the photoelectric absorption, in units of 10 21 cm -2 . 

t The photon index for the power-law model, and the value of p for the p-free disk model. 

§ An apparent innermost disk radius, assuming a distance of D = 3.4 Mpc and a disk inclination of i = 60°. 

Black hole mass estimated by assuming that the true inner radius i?; n coincides with 6i? g . 
Here the value of R ln is estimated by n n with correction factors of £ = 0.41 and n = 1.7. 
II The absorption-uncorrected flux, in units of 10~ 12 erg s — 'cm - 2 . 

# The bolometric luminosity, in units of 10 39 erg s" 1 . 

tt The 0.01-100 keV luminosity estimated from the best-fit model assuming isotropic emission. 
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Fig. 3. The 0.5-10 keV spectra obtained with the XMM-Newton MOS-1 on 2001 April 22 (left), and with the ASCA SIS (SIS0 
and SIS1) on 1999 April 6 (right). The best-fit MCD plus power-law models are superposed on the data. Bottom four panels show 
residuals between the data and the attempted models; from top to bottom, the MCD plus power-law, single power-law, single MCD, 
and the p-free disk models. 
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4. Discussion 

Using the ASCA data acquired on 1999 April 6 and the 
XMM-Newton data on 2001 April 22, we studied the X- 
ray spectra of M81 X-9. On both occasions, the source ex- 
hibited typical disk dominated spectra, with very similar 
sets of spectral parameters including the bolometric lumi- 
nosity of ~ 2 x 10 40 erg s _1 . Combined with the power-law 
type spectra observed from this object at somewhat (a fac- 
tor of ~ 2) lower luminosities (Ezoe, Iyomoto, Makishima 
2001; Miller, Fabian, Miller 2004), our results reinforce 
the reports by La Parola et al. (2001) and Sugiho (2003) 
that M81 X-9 makes occasional transitions between the 
disk dominant and power-law dominant states. 

In figure 4, we plot the present two measurements on 
the plane of Ldisk vs. T- m after Makishima et al. (2000). 
For comparison, we also show results on several high-Tj n 
ULXs, referring to Mizuno et al.(2001). Thus, the ULX 
data points, including those of M81 X-9, mostly fall on 
the apparent "super-Eddington" region, where Lboi sig- 
nificantly exceeds Le which is calculated assuming stan- 
dard accretion disks with R m = 6R g . This is just another 
representation of the "too high T; n " problem. The fig- 
ure also reveals the "variable r in " problem, because the 
ULX variations do not trace the dotted lines representing 
iboi oc TJ 4 , but instead, behave approximately as Lboi tx T? 
(or equivalcntly oc T^ 1 ). Thus, we consider that M81 
X-9 is subject to the same problems as described in § 1, 
although its variation is not very prominent on the dia- 
gram. 

In order to solve these problems with the standard- 
disk interpretation of disk dominant spectra of luminous 
ULXs, we have adopted the slim-disk scenario as a work- 
ing hypothesis, and analyzed the two data sets of M81 
X-9 using the p-free disk model. As a result, the XMM- 
Newton and ASCA spectra have both been described 
by this model more successfully than by the traditional 
MCD model. Furthermore, the derived temperature pro- 
files, p = 0.60 ±0.02 with XMM-Newton and p = 0.61+^ 
with ASCA, agree with each other, and are significantly 
smaller than the canonical value of p = 3/4 to be found 
with standard disks. In view of the theoretical prediction 
by Watarai et al. (2000) described in § 1, these results 
strongly suggest the presence of a slim disk in M81 X-9 
on these occasions. In particular, the good agreement be- 
tween the two different instruments significantly reduces 
the possibility that the deviation of p from 3/4 is an in- 
strumental artifact. 

In figure 4, results on the Galactic BHB, XTE ,11550 - 
564, are also presented for comparison. At luminosi- 
ties below ~ 0.4Le, the data points of this BHB thus 
trace tightly a constant-mass (or Lboi oc Tj 4 ) locus. At 
higher luminosities, however, the source starts deviat- 
ing from this relation (Kubota & Makishima 2004), ex- 
hibiting the "variable rj n " problem just like the ULXs. 
The same behavior has been observed from other BHBs 
as well, including 4U 1630 - 47 (Abe et al. 2005) and 
GRO J1655 -40 (Kubota, Makishima, Ebisawa 2001). 
Kubota & Makishima (2004) named this characteristic 



spectral state of BHBs apparently standard regime, be- 
cause the spectra in this state resemble those from stan- 
dard disks, except for significantly higher values of T[ n . 
They showed that the RXTE spectra of XTE J1550 - 564 
in this state can be described by the p-free model with 
p<0.75. 

Based on these considerations, we may presume that 
high- Tj n ULXs and Galactic BHBs in the apparently stan- 
dard regime have essentially the same spectral character- 
istics, which in turn are explained in terms of the forma- 
tion of a slim disk; the only difference is that the ULXs 
are by an order of magnitude more luminous. Then, the 
first problem (too high T- m ) can be solved at least qual- 
itatively, because a slim disk appears significantly hotter 
than a standard disk. This interpretation also solves the 
second problem (variable r; n ), because such a behavior is 
predicted by the slim disk theory (Watarai et al. 2000). 

Then, how about the third problem, too high transition 
luminosity between the disk dominant and the power-law 
dominant states? Luminosity differences associated with 
ULX transitions (a factor of ~ 2 in the particular case of 
M81 X-9) are generally too small to identify the tran- 
sitions with those between the low and high states of 
BHBs. In addition to the classical high-low transition, 
BHBs make transitions between the standard regime and 
the apparently standard regime at some critical luminosi- 
ties, which are typically several tens of percents of Le 
(Kubota et al. 2001; Kubota & Makishima 2004; Abe et 
al. 2005). Around this critical luminosity, yet another 
spectral state, called very high state (VHS; Miyamoto et 
al. 1991), sometimes appears. In this VHS, the spectrum 
becomes dominated by the power-law tail presumably due 
to Comptonization, even though the optically-thick disk 
emission underlies the spectrum Kubota et al. (2001; 
2004). Kubota, Done, Makishima (2002) argued that the 
power-law dominant state of ULXs is analogous to the 
VHS of BHBs, rather than to their classical low state. 
This idea was also applied successfully to other power- 
law type ULXs (e.g., Stobbart et al. 2006). Then, the 
ULX transitions between the disk dominant and power- 
law (plus soft excess) type states can be regarded as analo- 
gous to the BHB transitions between the apparently stan- 
dard state and the VHS. This solves the third problem, 
because the transition luminosity, ~ 10 40 erg s _1 , is now 
considered close to Le, rather than to several percent of 
L E . 

Finally, let us touch on the black hole mass in M81 X-9. 
If the observed spectra are securely confirmed as emerg- 
ing from the standard accretion disk, we can estimate the 
true inner radius, R m , from the measured, r ln , through 
the corrections employed in §3.1. We can then identify 
i?i n with the last stable orbit, to get an estimate on the 
black hole mass. However, now that the accretion disk in 
M81 X-9 has been confirmed to be significantly deviated 
from the standard disk, the correction from the apparent 
radius, r ln , to the true physical radius, R m , would become 
highly model dependent. Even if we somehow derived i?j n , 
we may not be justified in identifying it with 6R g in an 
extreme slim disk condition. We thus do not attempt to 
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Fig. 4. The calculated Lbol plotted against T; n , both ob- 
tained from the MCD + power-law model. The dotted and 
dot-dashed lines represent the relations of Lbol °c T.* and 
^bol T? n , respectively. The dashed lines are the relation 

of Lbol °c ^i^ 4 (constant Eddington ratio r/ = ibol/^Bdd)- 
All points of M81 X-9 refer to the present paper. Those of 
other ULXs, IC342 source 1, M81 X-6 and NGC 1313 Source 
B, are obtained with ASCA (Mizuno, Kubota, Makishima 
2001; Makishima et al. 2000). The data of the Galactic 
BHB, XTE J1550 - 564, are obtained with RXTE (Kubota 
& Makishima 2004). The disk inclinations of the ULXs are 
assumed as 60°. 
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Appendix. The analysis of SNR G21.5-09 

In order to examine the reliability of spectral and pho- 
tometric results on highly off-axis XMM-Newton sources, 
we analyzed archival XMM-Newton data of the supernova 
remnant G2 1.5-09, acquired in off-axis and on- axis point- 
ings. This source, one of Crab-like supernova remnants, 
has been observed many times with XMM-Newton for cal- 
ibrating the mirror vignetting (Lumb et al. 2004) . The ob- 
servations were performed on 2000 April 7 for 17 ks, 2000 
April 11 for 25 ks, and 2001 April 1 for 15 ks, with off-axis 
angle of 0.'08, 11 '12 and 13f02, respectively. We analyzed 
these data sets in the same way as for M81 X-9. Figure 5 
and Table 2 show the 1-10 keV spectra, and the param- 
eters obtained from single power-law fits to them. Our 
results on the on-axis data are consistent with those re- 
ported by Warwick, Bernard, Bocchino (2001). As shown 
in Figure 6, the off- axis spectral parameters are consistent 
with the on-axis ones and the absorbed 1-10 keV flux is 
reproduced to within 10%. 



estimate the black hole mass based on the radius param- 
eter obtained with the p-free fit. 

Instead of using the radius parameter, a much safer ar- 
gument on the black-hole mass can be made based on 
the luminosity. To keep the black hole mass to a normal 
stellar black hole of ~ 10-20 Mq, the observed luminos- 
ity is a factor 7-10 higher than its Eddington limit. Of 
course, such a highly super-Eddington luminosity could be 
possible from a theoretical viewpoint (e.g., Ohsuga et al. 
2005). However, it would be more natural to consider that 
the object is radiating at close to its Eddington luminos- 
ity, considering that highly super-Eddington luminosities 
have not been observed from Galactic/Magellanic BHBs 
even if the disk deviate from the standard disk supposedly. 
In this case, the black hole mass is estimated as ~ 100 Mq. 
This exceeds a typical upper end of black hole mass which 
can reasonably form from models of single stellar collapse, 
and may suggest some kinds of merging (Ebisuzaki et al. 
2001; Belczynski et al. 2004). 
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Table 2. Spectral fitting parameters for G21.5-0.9 (1-10 keV). 



offset 




rt 


Fx 1 


xVd.o.f. 


0.'08 

W.12 

13f02 


23.3 ±0.6 

22.4 ±0.6 
22.3 ±1.0 


1.86 ±0.04 
1.82 ±0.04 

79 +0.07 
x - ' z -0.06 


51.8 
56.7 
55.5 


446.5/434 
422.6/399 
298.8/323 



* Hydrogen column density for the photoelectric absorption, in units of 10 21 cm 2 . 
t Photon Index I\ 

t The uncorrected absorption flux in the range of 1-10 keV, in units of 10 — 12 erg s -1 cm~ 2 . 
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Fig. 6. The best-fit power-law model parameters of 
G21.5-09, plotted against the off-axis angle. From top to bot- 
tom, plotted are the f— fO keV flux (in f0 -11 erg s _1 cm -2 ), 
the hydrogen column density for the photoelectric absorption 
(in 10 22 cm — 2 ), and the photon index. The line drawn at 
12f5 indicates the off-axis distance of M81 X-9 owing the 
2001 April 22 observation. 



Fig. 5. MOS-1 spectra of G21.5 — 09 of acquired in; (a) an 
on-axis pointing, and off-axis pointings with offsets of (b) 
11'.12 and (c) 13'.02. The best fit power-law model is su- 
perposed on the data. 



